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1^ " A simple formula is derived for the Ricci scalar curvature of any smooth level set {i//(xo,a:i,. . . ,x„) = C} embedded in the 

^ , Euclidean space R""''-'^, in terms of the gradient Vi/a and the Laplacian Ai//. Some applications are given to the geometry 

' of low-dimensional p-harmonic functions and high-dimensional harmonic functions. 

S Introduction 

i Many problems in mathematical physics can be described by "inhomogeneous Laplace equations" in the Euclidean 
1*"^ I space Ayr(r) - V{r), which not only encompass the familiar examples of linear equations {e.g. the Poisson equation 
. in electrostatics, the Helmholtz equation in wave propagation, the Schrodinger equation in non-relativistic quantum 
mechanics) but also some important classes of nonlinear partial differential equations (including, but not limited to 
\^ the Landau-Ginzburg model for superconductivity, the Gross-Pitaevski equation for Bose-Einstein condensates and the 
Navier-Stokes equations for viscous fluids). 

In this series of works, we study the geometry of level sets associated with solutions to the "inhomogeneous Laplace 
equations" Axf/{r) - V(r), starting with the current article which revolves around a simple formula for the Ricci scalar 
curvature Sg{r) at any point on a smooth level set {y/(r) - C} embedded in the Euclidean space of arbitrary dimension: 



5«(r) = -Alog|Vi//(r)| -I- V ■ 



Aur(r) — ^ 

|Vi//(r)|2 



J> , The right-hand side of the formula above invokes the gradient Vi//(r) and the Laplacians Ai//(r), Alog|Vi//(r)| evaluated 
C^l ' in the Cartesian coordinate system of the ambient Euclidean space. Likewise, the divergence operator V- also refers to 
. the flat space formulation. By convention, we will restrict our attention to cases where the gradient is non-vanishing 
Vi//(r) ^ 0, so that the appearance of \Vy/ir)\^ in the denominator is meaningful. 

We call the proposed scalar curvature identity a "simple formula" because it only involves some familiar operations 
and quantities of physical interest. Especially, in many practical problems, the modulus of the gradient |Vi//(r)| provides 
a measure of the "field intensity", and Aiff{r) - Vir) is given as a prescribed function in the Euclidean space or as a 
I certain function(al) of y/(r). The simple formula in the preceding paragraph thus allows us to deduce some useful 
information regarding the curvature properties of the level sets and their relations to "field intensities", which we will 
• ^ describe in several papers in this series. 

In the current work, we open with generic formulations for the geometry of level sets (§ 91.111.21 1. and follow the 
^ proof of the aforementioned scalar curvature identity ( 91.31 1 with a brief discussion on its invariance properties ( 91. 4^ . 
- - - We will then show that the degenerate cases of the scalar curvature identity in low dimensions recover some familiar 
results in classical analysis ( 92. 1^ and are consistent with various worked examples in classical physics ( 92. 2^ . We 
close this proof-of-principle article with a simple low-dimensional application to p-harmonic functions (© and a high- 
dimensional application to harmonic functions (@. 



1 A Simple Scalar Curvature Formula for Level Sets 

1.1 Metric, Connection and Curvatures for Level Sets Embedded in Euclidean Spaces 

A natural way to describe the local geometric properties of a smooth level set {f{xo,xi,. . . ,x„) = C} as a hypersurface 
in IR"'"'^ is to introduce the line element for curves on the hypersurface as ds^ - gijdu' du^ , where a Latin index takes 
values from 1 to n, and every repeated index implies summation [Ref 1, §7.3, pp. 68-69]. Here, (u^,...,u") are local 
curvilinear coordinates on the surface (as opposed to the Cartesian coordinates (xo,xi,... ,Xn) e IR""^^), and (gij) is the 
(covariant) metric tensor. Superscript and subscript indices stand for contravariant and covariant components, respec- 
tively. Written explicitly, the components of the metric tensor takes the form gij - diV-djr, where 3; is a short-hand for 
d/du\ The contravariant metric tensor (g''-') is simply defined as the matrix inverse of(,gij). 
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The tangent vectors div at different points on the level set in question are connected by the Gauss formula [Ref [3, 
§30.4, p. 311]: didjr^T\.dkr+b ijfi, where n is the unit normal vector of the hypersurface, the connection coefficients 

are T^j - g^'^didjr-d/r and the coefficients of second fundamental form are bij - didjr-n. The connection coefficients 

can also be evaluated from the Christoffel formula [Ref 1, §29.3, p. 293]: T^^. = ^g'^^idigej + djgn - d(gij). 

The components of the Weingarten transform W - {b\) is defined by b\ - g-^^'bki and appears in the Weingarten 
formula: din - -b'^^djr, that is, dn - -WAr for infinitesimal changes tangent to the surface. The mean curvature is 
equal to ^ times the trace of the Weingarten transform: H :- ^ Tr(W) - ^g'^-'bij. 

The expression Rij - dfT^.j - djT^.g + T^^F^ - T^rj^ gives the iij)^^ component of Ricci tensor on the hypersurface 

[RefQ §37.4, p. 394], and the Ricci scalar curvature is defined by 5? = g''^Rij. 

For convenience, we shall borrow some notations and terminologies from physics: "force field" is defined via the 
negative gradient F{r) :- -Vi//(r); "field intensity" F :- \F\ is the modulus of the gradient; "field line" or "F-line" is the 
integral curve of the vector field F(r),re fS""""^. By convention, we orient the normal vector as n = F/F. 

To extend the discussion on an isolated level set (^-dimensional manifold) to a family of hypersurfaces that fill a 
region in the ambient {n + l)-dimensional Euclidean space, we can employ the natural curvilinear coordinate system 
Hu'^ ,u^ ,...,u"-) as follows: 

1. For every point r, the zeroth component of curvilinear coordinate coincides with the "level value": - y/ir), thus 
the coordinate is constant on any level set; 

2. Two distinct points r(i//,u^,...,u") and r(if/' ,u'^ ,. . . ,u'"-) can be linked by the integral curve (F-line) if and only if 
uj - u'j ,j - l,...,n, so the F-line acts as the curvilinear coordinate curve of Accordingly, F - -Vtf/ is equivalent to 
dor^-n/F^-F/F^. 

From the above curvilinear coordinate system constructed in a neighborhood surrounding a non-critical point r 
where F(r) - |Vi//(r)| ^ 0, one can define line element for curves in (n + l)-dimensional Euclidean space as ds^ - 
gfivdu^Au^ where a Greek index takes values 0,1,..., «. By definition, ^oo - F~^ - 1/g^^ and got - g^'' - 0, so the 
decomposition of the Euclidean metric is a direct sum of the "level value coordinate" and "level set metric". One may 
extend the definition of connection coefficients as d^d^r - T'l^^d^r, where the newly-arisen connection coefficients will 
be computed in the following proposition. 

Proposition 1.1 (Connection Coefficients) Suppose that a real-valued three-times continuously differentiable func- 
tion y/ e C^(S);IR) satisfies the "inhomogeneous Laplace equation" 

Ai//(r) = V{r), re^ 

in a certain domain J) of the Euclidean space IR""""^, and has non-vanishing gradients therein F(r) :- |Vi//(r)| ^ 0,r e 
J) c IR"'''^. We endow the level sets ofy/(r) with unit normal n - -Vi///|Vi//| and decompose the Euclidean metric ds^ - 
g^vAu^Au^ - F~^(diff)^ -\- gij du' du-' into a direct sum of the "level value coordinate" and "level set metric", then we have 
the following computations for connection coefficients involving the index 0: 

^0=-!^, r%4^ (1) 

ro._:^^ rj ^Ajm^^^JmdF^. (2) 
2 dy/ ' °° 2^ du"" F^^ du""' 

Tl^y'^dogoo^-^logF, (3) 
along with a modified version of the "harmonic coordinate condition": 

rO:=^'^-r°, = F2(r°o-r-o) = -v. (4) 

Proof To prove the three identities in Eq.[l] it would suffice to compare the equation df^dyr = T'l^yd^r with the Gauss 
and Weingarten formulae: 

-^-I—^r'^-^-Fb— —iF—]^b''-^ 
du'duJ du'' ''dip' duj\ d(p) J du''' 

The two identities in Eq.[2l as well as Eq.[3j follow from the Christoffel formula T^^ - ^g^^id^grjv + dvgfjTj - dr^gf^v)- 
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Juxtaposing the two expressions of in Eqs. [Hand 121 we can put down 

'■J 



dgij _ 2 d\ogdet(g,j) _ ijdgij 2 ij _ 2nH _ 2Tr(W) 

'W'F'j' 'W'F^ — 



On the other hand, the "inhomogeneous Laplace equation" Ai//(r) - -S/ -Fir) - V{r) prescribes the divergence of force 
field F(r), i.e. {F/y/g)do{F^) - V (hereafter g - detigij)), which can be combined with Eqs.[l]and[5]into 

as stated in Eq.[4l If y(r) = 0, and iff(r) is a harmonic function, then the formula - g^^T^y - hearkens back to the 
"harmonic coordinate condition" in general relativity. ■ 

Remark From the identity dogij - 2bijlF, we can also readily deduce dog^j - -2g'-^bpF. 

Another by-product of the foregoing argument is the following result that will be used later in this work: 

dF V . V 

Tr(W)+ — = i.e. (wV)logF = Tr(W)--. (6) 

difr F F 

When V - and n-2, the formula above is a standard exercise in electrostatics [Ref 2, Exercise 1.11]. □ 

1.2 Laplace Operator in Curvilinear Coordinates and Evolution of Mean Curvature 

The Laplace operator on the Euclidean space R""""^ can be presented in curvilinear coordinates as 

= g^'^'idf.d, - rf^yd;,) = ] 3^ (^^Vdet(^;,v)aJ . 

Vdet(g-^v) ^ ' 
Here, det(g'^v) - g/F^. Similarly, one can define the Laplace operator on the level set Z as 

Proposition 1.2 (Decomposition of Laplacian) The Laplace operator A can be rewritten as 

n 1 dF d d 

A^A^ + F^ — TT--^ r + V — , (7) 

dy/^ F^ du"^ duJ dy/' 

where 

Ai -^diig'jy^dj) = g'j{didj - r'ljdk) 
is the Laplace-Beltrami operator on the level set ofy/. Accordingly, we have the following formula 



A-. 



^ „ 1 r.29^'^ogF no ^1 n^^logF Tr(W)y 
AlogF = -FA^-^F^^^ tXo - -^^Az- _ ^ _. (g) 



Proof By definition, we have 



With the substitution of g^^ - F^ and the expressions for F^'^.,Foq,Foo,F° from Proposition ll.il we obtain the claimed 
result in Eq.[7l To derive Eq.[8l it would suffice to spell out the term F^Fqq using Eq.[6l ■ 

Remark By the relations Ar - and A^r - Tr(W)n, one may also use Eq.[7]to deduce an explicit expression for the 
second order derivative d^r/di/.^. □ 
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For a hypersurface embedded in Euclidean space, its Ricci tensor (Rij) and Ricci scalar curvature can be spelt 
out "extrinsically" in terms of the Weingarten transformation and the principal curvatures. Concretely speaking, in 
the principal curvature coordinate system where the Weingarten transformation W - (bj) is represented by a diagonal 
matrix with eigenvalues ki,...,kn {viz. the n principal curvatures of the hypersurface), one has 



iRij)i 



ijn^i,j<n 



'gll[ikl + --- + kn)kl-kl] Q 





gnn[ikl + --- + kn)kn-kljj 



■ iblbij - bkig''"'b,nj)l^ij^n- 



Here, to identify the two extreme ends of the equation above, we may recall that the geometric relation 

bMg''"'b,nj-b'lb,j^-R,j 

applies to any curvilinear coordinate system on any n-dimensional hypersurface embedded in IR""""^, as evident from the 
contraction Rij - R^^. of the Riemann curvature tensor R'^^j - bijb^ - biebK Consequently, the Ricci scalar curvature 

S& - g'-'Rij - b'^bj-bjbj can be explicitly given in terms of moments of principal curvatures - (Z"=i^^)^-ZJ5i=i^m - 
[Tr(W)]2-Tr(W2) 

In the proposition below, both the surface Laplacian Ax and the Ricci scalar curvature turn up in the evolution 
of the mean curvature along the "level value coordinate". 

Proposition 1.3 (Evolution of the Second Fundamental Form) We have the following identities 

^ = (6)6.,-a.a, + rf^.5,)i O) 

and 

aTr(W) 1 Tr(W2) ^ [Tr(W)]2-i% 

^-^-l-e; E^— = -^^iT ^ , (10) 

dif/ F F F F 

where is the Ricci scalar curvature of the level set in question. 

Proof From the identity doididjr) = diidodjr), we may deduce dor^j + TJjT% = diT'^j^ + T^j^^l^. This results in Eq. H 
upon substitution of the connection coefficients. Combining Tr(W) - g'-'bij and d^g'' - -2^**6-^ /F with Eq.[9l we obtain 

dTriW) 1 b^bi 

— T— - = -Azt; - -TT^, where 6^6{ = TriW^) = [Trimf-M. 
dip F F * 

This verifies Eq.[10l ■ 

Remark We originally discovered the scalar curvature identity mentioned in the introduction by an examination of 
the evolution equation for mean curvature (Eq.llOt. We leave it to the interested readers to recover such an "extrinsic 
proof". In the next subsection, we will only describe an "intrinsic approach" based on the definition of Ricci scalar 
curvature via the metric and connection on the level set. As will be explained elsewhere, the intrinsic method has 
better extendibility when the ambient space is not Euclidean, and/or the co-dimension of the submanifold is higher 
than one. □ 



1.3 Scalar Curvature Identity and Evolution of Ricci Curvature 

We are now ready to prove the formula 

Vi//(r) 



Sg(r) = -Alog|Vi/A(r)| + V ■ 



|Vi//(r) 



mentioned in the introduction. 



Proposition 1.4 (Scalar Curvature Identity) Let F - \ V^f \ and V - Axfr, then we have the following geometric iden- 
tity for the scalar curvature of the level sets for if/: 

dv 2Tr(W)y y2 „ a y y2 ivvw] 

AlogF + ^-- ^ + _.AlogF + ^-7r^-^-^.AlogF + ^-V.(^J.O. (11) 
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Proof We define the Ricci tensor for the local curvilinear coordinates (u^ -y/,u^,...,u"-)hy 

SO that g'^^R_^y - follows from the flatness of the Euclidean space K""''^. Meanwhile, by the direct sum decomposition 
of the metric gf^v, the following identity holds 

= g^'R,. - FHdxT',, - doT'oA + ro'orL - r^.r^,) + ^ + g^jidoi^ - djY% + rf/^o + 1^1^ - r°,r)o - r-r"^). (12) 

Using the modified version of the harmonic coordinate condition r" :- g^^T^^ - - V (Eq.|4ll, we may derive 

^ ^OO^At+S ^.^iJ^AO + ^O Om^~ ^OO^mT+S ij'^ mO + ^ '■^ 00 + ^ Om-' " "'^ ^ 00^ mi + S ^/j^mO '•^ 00 + ^ Om''- 

Combining this result with the identity '5or^'^. + FJ'^r^Q - diV'^Q + T^qF^., we may rewrite Eg. [12] as 

= FhdmTZ - doTZ + r- r^, - rs,ro\) + ^ + ^'^(rror", - + ^T/L - ^a^-o - ^To^'j- Vir'oo + C^) 

— JT^c 3 T^m _ p, T^m , -pm pT _ rm pA _ pO i _i_ rt2> _ „ij F"* _ Vr™ 

-r yUm^QQ "Oi Om ^ 00^ mi ^ OA^ Om ^ Om^ OO-" S ^ im^ jO *' ^ Om 



2iogF-3o(ro"5„-r«o) 



+ r '.Oml 00 + ^ 00^ m/ ^Om^OO-' + '^ ^ i 



m 

Om' 



where we have employed the relation F'^^'^g^Q^^ S^'^^^inJ-'jo = in the penultimate step. Exploiting the identities 
itemized in Proposition ll.il along with the relation - dm^ogs/g, we may reduce Ea.[T2linto 



= F' 



' 5 , dV 
logiT- 

01//^ 01// 



du" 



duJ Fj ' F^ 



dF 1 ^^aiog^^aiogF 



-.F^"logF-^ + F^ 



dV 



1 d 1-^ a 1^ igj"'dy^ d 1 



1 „, 



F du"" [ duJ Fj F y/g du"" duJ F 
dif/ 



dV_ 



which completes the proof ■ 

Remark In Eq. [Ill the expression dVld%r is understood as the derivative with respect to the curvilinear coordinate i/', 
thus is well-defined for any continuously differentiable inputs V e C^(Ti;U). □ 

Proposition 1.5 (Evolution of Ricci Curvature) We have the following formulae valid in arbitrary dimensions: 



1 mij-b^fRij 



k 



F 
1 



2b',^^^ + 2g'^b)\^,^k^-I\^^c^ 



1 - 1 



— a^irv^a.- 



ds) 

(14) 



where -.^ g'^Rij,^'^ := b^^g'^ - g'j b^j and G'^ ^ R'J - ^g'J is the (contravariant) Einstein tensor with R'j ^g'^Rl- 
Proof Using Eqs.[9land[T0lfrom Proposition |1.3l we may differentiate both sides of the identity = b^bij-bkig^'^bmj 
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to obtain 
dR 



A i 

F I 



(bMg'"^b 



mj 



-did,— + rf-dk- 



(be^g'Pbpk 



-didk- + Tf,d,- 



e'""b 



{bejg'Pbpk 



J " F Jk F 



S "mi 



1 

■F 



-b^b'' 

—^b-+b^ 



lblb,j-R,j 



■' F '■■I F 



2bk^g^'{bZbej-Rej) (b'^^.bik-R^k 



-dA-^Tl,d,-b^^ 



bZbjk-Rjk 



■bij^I.- + 



1 S&bij-b^^Rij 



+b^-^,^,-+T^^^^,-ub^^A^-^l,^,-ub^^,^,--^^^^^^ 



thus proving Eq.[T3] (Here, we have used the relation b^b^j - bkig b^j - bjb^i, which also implies the symmetry 

b\Rej^b'.Rei.) 

Carrjang this further, we have 



m_ 2gi>'biRij ..dRij 
dy/ F dy/ 



2biR] 



.2b',^^^ + 2g^^b)\^,^UJ-T\,^,^\, 



as claimed in Eq.[T4l Using covariant derivatives indicated by semicolons, we may rewrite the right-hand side of Ea.[14l 
as 



2biR] 



2biR] 



-2\fi'^dij\^+2(fi"'),kdi^^- 



2biR) 



—di\p'''^du^\+2(p"'),kdii; 



but 



= b',.^g'^ -g'^b].^ = b',.^g'''-g'^bij = 0. (15) 
This leads to a succinct proof of Ea. [T4Fl based on the fact that covariant derivatives of the metric vanishes (g''*)-^ := 

. . :\b^.-Yf..b) _ 



dfgik + g'"'T^ f + g''"'T' . = 0, and the symmetry {dib^. + r*,6^ - V''. .b) = 6*. . which follows from the Codazzi- 



Mainardi equation dib^j + T^^b^j = djb] + r*^6 



1.4 Diffeomorphic In variance of the Curvature Formulae 

Locally, a level set of yi can be identified with a level set of f oy, where the smooth function f : IR — ' R has non- 
vanishing derivative f ^ 0. It is thus natural to ask if the formulae for mean and scalar curvatures (Eqs. [6] and llll l 
would remain intact under such function compositions. This question is answered in the positive by the following 
proposition. 

Proposition 1.6 (Diffeomorphic Invariance) For every f £ C^(IR;IR) satisfying f ^ 0, Eq. |6]is diffeomorphically in- 
variant in the sense of 



d/-/di//^^^^^^ aiv(/-o,//)l _ Mfoyj) 



\AflAyj\ 



d(foy/) \V(foy/)\ 



(16) 



where Tr(W) is evaluated with respect to the normal n - -Vy//\Vy/\, and the diffeomorphic invariance of Eq. [77] is 
manifested in the following way: 



Alog|V(foi//)| + Sg- V- 



Mfoy/) 



v(/'°y) 

|V(/-oi//)|2 



:0. 



(17) 



Proof By the chain rule, we have 



V(/-oi//)=l/^Vi//, Aifof)^^Af+^\Vy/\^, 
ay/ ay/ di//^ 
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thus we may compute 

A(/-oy) a|V(/-oy)| 

IVCZ-oi//)! d(foyj) Idf/di/'l IVi/'l ' Idf/di/'l 
which confirms Ea.ll6lthrough Eq.[6l Next, we evaluate 



df/dy/ dyj 



dy/ 



dfldyj I Ay/ d\Vy/\ \ 
\df/dy/\ iivi/^~ dy/ )' 



Aifoy/) 



Vifoy/) 
|V(/-oi//)|2 



-Alog|V(/-ov/)| = V- 



df d^f o 
ay/ dyf^ 



d^f/dyr^ o 



dv 



df/di/' 



-^Vy/ 
dy/ 

Vy/ 



Alog 



dy/ ^ 



\Vy/V 



-Alog 



di// 



V- Ai// 



■ AloglVi/'l + V ■ Vy - Alog 



and quote the decomposition of Laplacian (Eq.O to deduce 



^ d^fldy/^^ 

V- — ^Vw 

I d/Zdi/A ^ 



-Alog 



d/- 



di/' 



■.\Vy/V 



dy/' 



■log 



di/'' 



■log 



dy/ 
df 



dy/ 



+ (Ai//)— log 
dy/ 

d 

+ (Ai/')— log 
di// 



d/ 



dy/ 
df 



Alog 



d/'/di// 



d/- 



df 



dy/ 



dy/ 



dy/ 
32 



\^y/r-^ + {Ay/)— 



d 

' dyr 



log 



df 



dy/ 



which completes the verification of Eq.flTl 



2 Scalar Curvature Formulae in Low Dimensions 

2.1 Examples from Classical Analysis 

The Ricci scalar curvature identity (Eg.lllll in low dimensions can also be independently verified with brute force. 

Example 2.1 (Ricci Scalar Curvature in Low Dimensions) For any real-valued, three-times continuously differ- 
entiable function yi e C'^{Q.\K) defined in a two-dimensional domain Q c K^, one has the exact identity 



(dly/ + d^yy/)dxy/ 
(d^y/f + (dyy/f 



{dlyj + d^yy/)dyy/ 
(d:^y/f + (dyy/f 



wherever |Vi//| 7^ 0, 



(18) 



which degenerates into the obvious relation d'^\og\dxy/\ - dx(d^y//dxy/) when there is no dependence along the y-axis. 
For a surface embedded in depicted by the implicit function y/(x,y,z) - 0, the Gaussian curvature K - Skl2 satisfies 



2K^dx 



idly/ -\- d^y/ -\- d^y/)dxy/ 



+ (3^ 



(3^1//)2 + (3y1//)2 + ((3z1//)2 

(dl + dl + d^) log J{dxy/)^ + (dyy/)^ + {d,y/)^. 



(dly/ -\- d^y/ -\- d^y/)dyy/ 
(dxy)^ + (dyy/)^ + (d^y/)^ 



+ (3, 



{d^y/ -\- d^y/ -\- d^y/)dzy/ 



{dxy/}^ + (dyy/)^ + (d^y/)^ 



In particular, for the Monge form z — f(x,y), one has 



[l + ((3^f)2 + (5„/-)2]2 



dlf dxdyf 

dydxf dlf 



(19) 



(20) 



Proof Writing y/ix,y) as a bivariate function oi ( - x -i- iy,( - x - iy, then we have d^y/ - (d^y/ - idyy/)l2 and d-^y/ ■ 
{d^y/ + idyy/)l2, and both sides of Eg.fTsleaual 23^5^1og(3^i//5^). 
The right-hand side of Eq.[T9]is equivalent to 



[(dxy/1^ + (dyy/f + {^^y/nH^^y/f 



X det 



d^y/(d'iy/d^y/ - 2dxy/dxd^y/) + (d^y/fd^y/ 



d^y/i-d^y/dyd^y/ + d^dyy/d^y/- d^d^y/dyy/) + dxy/dyy/dly/\ 



[dzy/i-dxy/dydzy/ -\- d^dyy/d^y/ - d^d^y/dyy/) + d^y/dyy/d'^y/ d^y/idyy/d^y/ - 2dyy/dydzy/) + {dyy/) d'^y/ 

which is a standard formula in differential geometry. Taking y/{x,y,z) -z- f(x,y), we recover Ea.[20l 
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Remark Despite being a vacuous truth, the differential identity Eq. [18] (when appropriately rewritten) can serve a 
good purpose for the geometric analysis of two-dimensional problems. See, for example, the applications in ^ □ 



Example 2.2 (Planar Harmonic Function) Let (p(x,y) be a harmonic function in a certain domain in R , then 



, 

^ + ^)log|V(p(x,3')|=0. 



Proof From complex analysis, one can rewrite \V(p(x,y)\ - \ f'(.0\ where ( - x + iy and f(0 is constructed from (p(x,y) 
and its harmonic conjugate. As log\V(p{x,y)\ - Re\ogf'(0 is the real part of a complex analytic function, it must be a 
two-dimensional harmonic function. ■ 

Example 2.3 (Conical Harmonic Function) In spherical coordinates r — (rsin0cos0,r sin sine/), r cos </>), the func- 
tion (p(r) - /i(tan |,0) is harmonic so long as h{p,(p) defines a two-dimensional harmonic function in polar coordinates: 
((5p -I- p~^dp + p~^d^)h(p,<p) - 0. The level set of <p(r) is a conical surface (necessarily with vanishing scalar curvature) 
excluding the apex r-Q. The following identity is satisfied for \r\ ^ Q: 



Alog|V^(r)| = -5- 



— ir'^— + 
dr \ drj sin (30 



d ( d\ 1 

sm0 — H T 

de] sin2 0#2 



log|V(p(r)| = 0. 



Proof With the substitution p - tan |, we obtain 



,j ^2 2.^.^2 (d0)^ + sin^0(d0)^ Id Id' 

{Apr + p'^(A(pr^ — — — and — -H- 



4 cos* 



dp' p dp p' d(p' 



■ 4 cos 



1 ^ ■ 1 
^— — sine— -H 
sine de[ dO sm 



^ed(p' 



Therefore, we have 



3/i(tan|,0)' 

de 



sin^e 



5/i(tan|,(/))' 



d(l) 



(dh(p,(l))Y^ 1 {dh(p,(p) 



dp 



p2 \ d(p 



and 



Alog|V(p(r)| = -5- 



d o d] 1 d ^d 

— k — + sine — 

dr[ dr sinO dO \ dO 



log 



2rcos2 ' 



1 [ d' ld_ J_d^], IdHp^Y 1 l dh{p,cl)) 
4r2cos4|i(3p2 + p(3p^p2(302j "Syi dp j ^ pA dcl> 



d^ Id 1 (52 

-I- h 



jl dh(p,(p) \^ ^ 1 l dh{p,(p ))^ 

4/-2 C0S4 I I (5p2 ■ pap ■ p2(5(/)2 J °^V' dp j p2[ (50 



Now, quoting the result in Example l2.2l we see that Alog|V(p(r)| = 0. ■ 

2.2 Examples from Electrostatics 

Many worked examples in electrostatics provide analytic expressions for the geometry of equipotential surfaces 
(level sets of the electrostatic potential (p{r), which is a harmonic function) and the corresponding spatial distributions 
of the electric field intensity E{r) - \V(p{r)\. Several specific cases listed below thus serve as direct verifications of the 
relation Alog£(r) -i- 2K(r) - 0, where K{r) - 5?(r)/2 is the Gaussian curvature of the equipotential surface. 

Example 2.4 (Isolated Ellipsoidal Conductor) For an isolated conducting ellipsoid surface 



x' y2 2^2 

^ + 7^ + -5- = l {a>b>c) 



(21) 



with total charge Q, the surface charge is distributed as [Ref\^ p. 22] 



crir)-- 



4nabc 



x2 y2 z' 



4 64 c4 4n V a262c2 



Q 4 K{r) 
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thus on the conductor surface (with r{x,y,z) fulfilling Eq. [2i] . the electrostatic field intensity satisfies \ogE - log + 
const and the following identity is satisfied: 

Proof Following the convention in Ref. 13], we introduce the curvilinear coordinates r - ri(,T],() > -c^ > ry > -b^ > 
( > -a^) as the three roots of the following cubic equation in u: 



a^ + u + u c^ + u 

From the argument in [Ref Is^, §4], the electrostatic potential is given b50 

-+00 



2J( 



k y^{s + a^){s + b^){s + c^) 
and field intensity is distributed as 

logi;«, r?, C) = - log ^/(.^-7lM-0 + const, 

and the equipotential surfaces (the ^ - const surfaces), being ellipsoids confocal with the one shown in Eq. |5T1 has 
Gaussian curvature 

(a^ + 0(b^ + 0(c^ + 

^'^^'^'O = 77 ,2.f ;■^2 • 

i(-T]Y{(-0 

With the Laplace operator [Ref 0, p. 19]: 



4 



ii-m-oiri-o 



(,_0|E,-j H(-0[R,^^] -^^-^Y^e-, 



Ru^V(u + a^)(u + b'^)iu + c2), u = ^,T/,C 

one can verify that 



(C-??)2(e-02 



indeed holds. 



Remark The ellipsoidal coordinate system employed in the proof becomes ill-defined for the spheroid case (where 
either a - b > c or a > b - c), but the relation E oc \/K and 2K + Alog£ - remains valid, as the readers may check on 
their own. 

In the fully degenerate case a-b-c, the equipotential surfaces form a family of concentric spheres. On the surface 
of radius r, the Gaussian curvature K - 1/r^, and the spatial distribution log£(7-,0,0) - -21ogr + const. It is trivially 
true that 

AlogE- ——ir^^^^^]--—--2K 
r^ dr\ dr j r^ 

For a "point charge" in higher dimensional space IR""'"^, we have a power law log£(r) - -nlog|r| + const instead of the 
Coulomb inverse square law. Then we may perform the following computations for a hypersphere {x £ IR""'"^||a:| - r} of 
radius r: 

Al E ^^^^^^ n dlogE n(n-l) 
Qj.2 r dr ' 

which is consistent with the fact that a hypersphere of radius r in IR""""^ has scalar curvature - n(n - l)lr^ . □ 
Below are two more examples concerning conductors at electrostatic equilibrium with external sources. 



^We have converted the statements in [Ref. 3, §4] to SI units, wherever apphcable. 
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Example 2.5 (Rectangular Box Conductor Interacting with Electric Charges) Consider a grounded rectangu- 
lar box conductor at electrostatic equilibrium with a distribution of electric charges (with density g{r)) that vanishes in 
a neighborhood of the box boundary: 

"^^(P^r) = r £ no := (0,a) x (0,6) x (0,c),supp£) d Qq; 

[(p{r)-0, redQo- 

For a boundary point r e SHo that is not situated on the edges or comers of the rectangular box, we have Alog|V(p(r)| - 0. 

Proof The electrostatic potential ^(r) can be represented as eo(pir) - J/Jq^ p(r')G£)(r,r')d^ r' where the Dirichlet Green 
function GD(r,r') is given by 

g oo oo ^ sin ^ sin ^ sin ^ sin ^ sin ^ sin ^ 

a'' 

Now, at a boundary point r e (0,a) x (0,6) x {c} on the top side of the rectangular box, we have d^(p(r) - d^(p{r) - d^(p{r) - 
dxdy(p{r) - judging from the properties of GD(r,r'). In order to verify that Alog(£^) - 0, it is sufficient to check 

E^A(E^) = 2E^Tr{[VV(p(r)f} = 4[d,(p(r)f {[dxd,(p(r)f + [dyd,(p(r)f} 

according to the Bochner-Weitzenbock formula in Euclidean space, and 

|V(£2)|2 = {aj|a,(p(r)|2]}2 + {dy[\d,ip(r)f]f = 4[d,(p(r)f {[dxdMr)f + [dydMr)f} 

for r e (0,a) x (0,6) x {c}, the identity Alog(£2) = is verified. Similarly, we can double-check Alog|V^(r)| = on the 
other five faces of the rectangular box. ■ 

Example 2.6 (Ellipsoid Conductor in Uniform External Field) Consider a grounded conducting ellipsoid surface 
with axes a> b> c (defined by Eq. \2D , placed in uniform external field EqCx- The electrostatic potential [RefJs, p. 23] 



i: 



ds 



, , „ (^ + a2)(r? + 62)(C + c2) ^^ Jj {s + a^W{s + a^)(s + b^){s + c^) 



(s + a2 



i2)V(s + a2)(s + 62)(s + c2) 

satisfies 



^a^b^c^ 

on the conductor surface where ^-Q. 



Alog|V(/)(r)| + 2— — - = 



Proof Upon the substitution F(0) - 1, we have the following expression for ^ = 0: 



2 I a262c2 



Alog-E(r)= — " " -2a^6^c^ 



F"m 
F'(0) 



' o 2^2 2-P""(0) 2,u2 2,F"m (a^ - 3c^){b^ + C^) 2 o ^'^^ 

+ 2a^b'^c'^ , + a^(b'^ + C^)—— 5 +a^- 3 — 5- 

F'(0) F'(0) 2c2 a2 



With the relations 

F"m _ 3b^c^ + a^ib^ + c^) F"'m _ 15b^c^ + 6a^b^c^{b^ + c^) + aH3b^ + 2b^c^ + 3c^) 
F'(0) ~ 2a^b^c^ ' F'(0) ~ 4:a^b^c^ 

we obtain 

a%h^ 

AlogE(r)^-2^-^, 

which is the claimed identity. 
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3 A Simple Application to p -Harmonic Functions in Low Dimensions 

We consider a class of nonlinear partial differential equations 

V-(|Vv/|^-2Vi//) = 0, p>l. (22) 

A solution xf/ e W^'P(0.) (2] to the above nonlinear equation minimizes energy functional /q |Vi/a|^ d'^r,d > 1 (which is 
the LP analog of the Dirichlet form /q |Vi//|^d'^r) with homogeneous Dirichlet boundary conditions on each connected 
component of dD.. Functions satisfying Eq.|22]are referred to as p-harmonic functions tsll. By convention, the attention 
is restricted to the scenario with nowhere vanishing gradients S/xf/{r) ^ 0,r e Q, so that one may recast Eq. |22]into the 
form of "inhomogeneous Laplace equation" 

dF 

= Ai//+(p-2)Vi//-Vlog^' = Ai//+(p-2)F — , where F -.^ Wy/l. (23) 

aifr 

Now that V = {2- p)FdF/dy/, we immediately arrive at the following geometric identities, according to Eas.[6land[TTl 

Tr(W) = (l-p)^, (24) 
ay/ 

c,d^\ogF oldFf 1 ^a^log^' IdF]^ 

^.-Mo,FH2-pW^^H2-pf[-] .FA,-.a-pW^^Hl-pX2-p,[-] 

1 ^ ^.^.n-nW^f 



.^A.-.(l-p)ir— .(l-p)^^-J . (25) 

Here, Eq. m is consistent with Eq. 4.3 in Ref 15]. From the geometric inequality id - 2)[Tr(W)]^ -id- 1)S& ^ id - 
l)^var[spec(W)] > 0, we obtain 

For p ^ 2, we may rewrite Ea.l25las 

1 l-o d^iF^'P) 1 d Tr(W) 

F 2-p dy/^ F dy/ Fp-^ 

and the rightmost expression in the equation above remains valid for p-2. 

Now, we may test the validity and usefulness of the formulae above in the following proposition by recovering the 
d = 2 case of Theorem 1.2 in Ref dl. 

Proposition 3.1 (Extremal Properties of Planar p-Harmonic Functions) Let y/ir),r e D.cU'^ be a p-harmonic 
function, with p>l. Then its level sets satisfy 



d S& _ 1 ^ 
'dyj^P ~ FP-^y^ ' 



TriW)\ 2 



Tr(W) 



(28) 



where ibj^) is the matrix representation of the Weingarten transformation W. 
In particular, for d-2, we have the equality for -F^~PTriW) - \Vyf\^~PK 



\ K K ip-2f l dlogF f 



+ F'^\-—^+2 _ ' \^^\ =0, (29) 



wherever ViF^'Px) - 0, so that a non-negative-valued |Vi//|^~^k attains its minimum at the boundary of a domain. (Here 
in Eq. d/ds is the tangential derivative with respect to arc length parameter of the level set, and Ax = d^/ds^ is the 
Laplace-Beltrami operator on the level set.) 
Furthermore, for d-2, we have the equality 

a( Ai I f 2 ] K _ a(a + 2-p) k 2(p - 2)^ + a:(a + 3p - 6) / aiogF \^ 

p-1 df^j FP-'^+" ~~ FP-^ (p-l)2 FP-i (p-1) i ds 1 ' 
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wherever V(F^ p ")c) = 0. Thus, if either of the following 



a£(-cx),p-2]u[0,2-p]u[2(2-p),+cx)), l<p<2 
a e (-cx),2(2 - p)] u [2 - p, 0] U [p - 2, +oo), p>2 



(31) 



happens, then a non-negative-valued |Vi//|^ p "k attains its minimum at the boundary of a domain. For 3/2 < p < 3, the 
minimum of a non-negative function K{r),r eDu dD is attained at the boundary dD. 

Proof We divide both sides of Eq.[27]by FP and differentiate, to arrive at 



d S& d i 1 1\ Tr(W) 



dy/FP dy/\FP-'^ F] dy/^ Fp-'^ Fp-'^^ dyj\ 



Tr(W) 



F] dy/^ FP-'^ 



where the p-harmonicity has been exploited in the form of d^{FP ^ ^) - 0. Now, we simplify 



with the relations dy,g'j = -2^'* and d^iF'PIp) = F'-P-^ Tr(W)/(p - 1), to obtain 



dy/ [ 



FP) p - 1 



giJFP-^^d 



Tr(W) \ 



which proves Eg. 1281 a formula not found in Ref [5]. 

For d = 2, we have S^ = 0, and g'''bl = g'-' Tr(W) = -g'-^K. This brings us from Ea.[28lto 



= 



FP-^Vg 
1 

F^VS 



'TriW)FP-^^dj — 
r P 



Tr(W)^ oP-2f ;,Tr(W) _^ 1 



32 Tr(W) 
' 'dyJ^ FP-'^ 
fl2 Tr(W) 



1 Ai(Fi-'PTr(W)) p-3 dlogF d Tr(W) p-2Tr(W) d Tr(W) 
(p-l)F2 ds FP-'^ p-1 F dy/ FP-'^ 



F^ p-1 

p-2 d 
- 2 - 

(p-i)FP-i as 



Tr(W) 



FP- 



d 1 ^ Tr(W) 



dsF dy/^ FP-'^ 



1 AzCF^-PTrCW)) Tr(W) 

+ 



ir2 



p-1 



dy/^ FP-'^ 



3p-7 aiogJ' a Tr(W)^2^-2Tr(W) 3 Tr(W) 



(p-Dii^^ 5s ds FP-'^ p-1 F dyr FP-'^ 



\ I \ VV I \ U — ^ J t U llthH' ^ ^ 

+ 2- 



Tr(W) (p 


-2)2 


/aiogF 


FP-^ (p- 


1)F2 


i ds 



(32) 



where s is the arc length of the level set, and we have evaluated Aiil/F) from Eq. [27] (with = 0) in the penultimate 
step. 

Now, for any real number a, at points where V[F^~p~" Tr{W)] - 0, we may transform Ea.l32linto 



= 



1 AziF'^~P~"T:riW)Fn Tr(W)^'" 
h 



F2 



p-1 



gyj2 pp-l+a 



3p-7 aiogF d Tr(W)F" ^^p-2Tr(W) d Tr(W)F'' 



(p-l)F2 ds ds FP-'^^" p-1 F dyr FP-'^+" 



+ 2 



Tr(W) (p-2)2 I dXogF V 
FP-'^ (p-l)F2 i ds ] 

F Ai(i^ ^ Tr(W))- Az^ + a(a+l)-^|-g^J ^2 rp^.^^^ rp^.^^^ ^2(^a) 

h i*^ — h 



F2(p-1) 

a(3p-7)Tr(W) MlogF \2 ^ ^p-2[Tr(W)]2 5(^'«) 



5i//2 



-l+a 



p-1 F''+« dy/ 



+ 2 



Tr(W) (p 


-2)2 


/31ogF 


FP-i (p- 


1)F2 


i ds 



pp-l+a 
.2 



3l//2 



Aj;(iri-P-«Tr(W)) ^„ ^2 Tr(W) 



F2-«(p-l) 



(3i//2 FP- 



F2(p-1) 



a[Tr(W)]2 a(F«) a[Tr(W)]2 5(F2(«-i'+P) 



pp+a-2 p2(p-2+a) 



dy/ 
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Tr(W) 2ip - 2)2 + a(a + 1) + a(3p - 7) (dlogFY 



FP- 



A^iF'^-P-^TriW)) ^„ Tr(W) 
+ F 



F^-Hp-1) 



FP-'^+" 



[ ds ) 

[Tr(W)]^ - 2a{a -l)-ap + 2a{p - 2) 



Tr(W) 2(p - 2r + aCct + 1) + a(3p - 7) / aiog^' \ 



FP- 



{p-l)F^ 



\ ds 



as claimed in Eq.|30l Solving the inequalities a(a + 2-p) > and 2(p-2)^ + a(a + 3p-6) > 0, we obtain Eq.[3Tl Setting 
a = 1-p in Eq.[31j we reach the range 3/2 < p < 3. The conclusions drawn from Eqs. l29ll3T] agree with Ref. tSJ- • 

4 Application to Einstein-Hilbert Actions of High-Dimensional Equipo- 
tential Surfaces 



During the study of a physical problem that will be described elsewhere, we encountered the following quantity 

'.dmCLip) 



f 



which is the Einstein-Hilbert action of a curved equipotential surface (i.e. level set of a harmonic function with level 
value (/?). For dim(Zg,) - 1, the Einstein-Hilbert action vanishes identically. For dim(Z^) = 2, the Einstein-Hilbert action 
is equal to 4n times the Euler-Poincare characteristic X^^qiX so it remains a local constant of (p so long as the field lines 
are not broken. 

In the following proposition, we will present a convexity result for the less trivial cases where dim(Z^) > 3. To 
emphasize that we are handling harmonic functions, we will replace the earlier notation F - |Vi//| with E - {Vcpl. Here, 
(p(r) can be regarded as a high-dimensional "electrostatic potential", and E(r) = \V(p(r)\ the corresponding "electric field 
intensity". Still, it will be tacitly assumed that E - {Vcpl is non-vanishing. 



Proposition 4.1 (Conditional Convexity of Einstein-Hilbert Action) The Einstein-Hilbert action 
satisfies 



lbZG'Jbij-R'jRij + ^^%^ 



dmCLJ-f Gikdi—dk — dmCZa,). 

la, ■O' E 



6 £2 
Accordingly, for n>3, we have s^"{(p) > i/'Z^ has non-negative sectional curvature. 
Proof It follows directly from Ea. ll4Fl that 



(33) 



G'Jbi 



■dm(Z^) = 2 
E 



Gijip^j-bW') 



E 



dmCLip) - 2 



E 



dm(Z, 



E 



■dm(Z^), 



where Gij - Rij - \gij^ is the covariant Einstein tensor. Before computing the second-order derivative s^"((p), we need 

:fik 



±d(P^ b_^ d_ , 

^ dip E ^dcp^^"'^ ^ ^ 



m oik 



bZP 



E 



bib';\ 2bZgijb] 2g''-bib] 2g'"^bib 



bZib'^g'^-g'jb)) 



-,ik 



E 

1 b^ b"^ 



E 



E 



Ibjbme 



E 



Ej-je 



Ax — + 
E E 



2b'^gijb) ^ 2g''^blb) X J2g^'^bib) g''bjbi^ 



E 



E 



urn iji^k 



E 



r)„imTfk a'^nut ui 



E 



E 



■^'>^E-E 



1 ^\ 3R 



•>ik 



E '^[e 



:ik 



^ ' ^E 2E 



3G 



ik 



E '^[e 



:ik 



(34) 
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Thus, we may evaluate 



diph^ E ^ h,^ E ^ d(p ^ dip E 

(n \ r S& 1 r R'jRij 5-n C S?^ 



d Gi 



where we have used the squared modulus of the Einstein tensor G'^-'Gij - R'^-'Rij + (f - in the last step. From 
Eqs.[T3]and ll4Fl it is not hard to see that 



d Rij _2Gik 1 1 



'2 6^ 



r,ij " ^'■IJ _ ^^ik ' _(a Rik r-. ^ 



-Dbl 



so we have 



dcp E £2 



2 d(p E £2 



— <f> — — ^dm(Z„)=(l) — r- — r-^^ — dTO(2™)+(i) /3'-'- — -fdm(Z^) 

dtpJz^ ^ E ^ dip ^ dip E 



d.(Z,) - 3/ ^ d.(Z,) . ^ £ ^ d.(Z,) 



+ 



as a consequence. Meanwhile, 



2G'jb 



'ij 



can be combined into the previous results to yield 

n\ d 



dm(Z^) + £ ^[-Az^ + ^]dm(Z, 



^ ^ d^i — - 2) d^i — <i-(^.) 



R'JR^ 



After we apply the contracted Bianchi identity (Gik)'^ - 0, Ea.[33]emerges 
In the principal curvature coordinates, the Einstein tensor behaves as 



'-^11 H klkr 



(Gij)i<iJ^n — 











-gnn ^^^n 



which is clearly a negative semidefinite matrix if has non-negative sectional curvature. In this case, we then have 

s^"iip) r lbZG'^bij-wmij + \<5t^ 



6 -Jx, 

However, we may use principal curvatures to rewrite 

hlG'jbij-R'jR.j + ^St^^ X kikjk] + S Y. kikjkekn,, 

" i<! i<i<(<m 



(35) 



(36) 
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where the second summand is left out for n -3. Such an identity as Eq. [36]is elementary though probably not self- 
evident. By inspection, the left-hand side of Eq. [36]is a homogeneous quartic polynomial in the variables (ki,...,kn), 
where each variable cannot appear in a power higher than two. Terms like k'^gk^ (where ( <m) will not count towards 
the net sum, because -R^-'Rij contributes -2k^k^ while ^5?^ yields 2k^k^ in their respective developments. Therefore, 
by symmetry, the left-hand side of Eq. [36] can be put into the form 

a kikjk^f + b kikjkekm, 

i<j i<j<e<m 

for some constants a and b. To determine that a-1, simply set ki-k2-k^-l and ki - ••■ - To find out that 

6 = 8, one may consider the special case where k\-----kn-'^- 

To summarize, Eos. [35l and [36l explicitly demonstrate the non-negativity of ^"{(p) for equipotential surfaces with 
non-negative sectional curvature. ■ 

Remark In the proof of the proposition above, we have used the following alternative forms of Eqs.[6land[T0l 

. dE dTriW) 1 [Tr(W)]2-Sg 

Tr(W)+ — = 0, — ^ = -Az--^-^^^ . 

d(p dip EE 

Here, the first formula is sensitive to the harmonic condition A(p(r) - 0, so the conclusion of this proposition does 
not apply to level sets of arbitrary functions. The second formula displayed above is morally equivalent to the scalar 
curvature identity 5? + Alog£ - 0, as can be seen from the one-line proof: 

1 od'^^osE 1 , d Tr(W) 1 aTr(W) , . d 1 

^ + Alog£ = ^-£Az-+£2 %-^S&-E^Y,--E^ -—^Si-E^z--E —-E^TriW) , 

E dcp'^ E d<p E E d(p dcp E 

which hinges on a degenerate version of Eq.[8l □ 
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